An analytical treatment of in-plane magnetotransport in the Falicov-Sievert model 
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We derive an analytical expression which allows efficient computation of the effect of all the Fermi 
surface trajectories induced by a combination of Bragg scattering and magnetic breakdown on the 
in-plane components of the resistivity tensor. The particular network of coupled orbits which we 
consider was first formulated by Falicov and Sievert, who studied the problem numerically. Our 
approach, based upon a method used previously to derive an analytical solution for interlayer trans- 
port, allows us to show that the conductivity tensor can be written as a sum of a matrix representing 
the effect of total magnetic breakdown and one representing a combination of complex electronic 
trajectories, and we find a compact expression for the in-plane components of the resistivity tensor 
that can be evaluated straightforwardly. 
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I. INTRODUCTION 

Magnetoresistance has long been used as a tool to 
study the Fermi surface of metals^. The observation 
of quantum oscillations provides a measure of the cross- 
sectional area of closed pockets (due to the quantization 
of the Landau levels). The background magnetoresis- 
tance also contains information about all the electronic 
orbits induced by the magnetic field. For example, a 
transverse magnetoresistance saturating in high held can 
be evidence for the presence of closed orbits, while a mag- 
netoresistance that increases quadratically with field can 
indicate open orbits. A number of effects can compli- 
cate this picture and one such is magnetic breakdown, a 
concept first introduced to explain a giant orbit observed 
in magnesium 3 . Because of the periodic potential in a 
crystal, small gaps in the dispersion sometimes open up 
at the Brillouin zone edge, splitting the Fermi surface 
into distinct sections. In low magnetic fields, electrons 
traverse the Fermi surface under the action of the mag- 
netic field and Bragg scatter at the Brillouin zone edge. 
In a sufficiently strong magnetic field it is possible for 
an electron to tunnel from one segment of the Fermi sur- 
face to another— This is because in high magnetic field, 
mixing of states between different pieces of the Fermi sur- 
face can increase the likelihood of magnetic breakdown, 
so that the resulting orbit more closely resembles that 
which would be obtained in the absence of the periodic 
potential^. 

The presence of magnetic breakdown presents a chal- 
lenge for the simulations of magnetoresistance. This is 
because at the points on the Brillouin zone edge in- 
tersected by the Fermi surface (points which we call 
magnetic breakdown junctions) there is a probability p 
for magnetic breakdown to occur and this is given by 
p = exp(—Bo/B), where Bo is the characteristic break- 
down field; the problem is that under usual experimental 
conditions p is between and 1 . Thus there is a hierarchy 
of complex trajectories that must be summed to account 



for all possible contributions to the conductivity in which 
magnetic breakdown either does or does not occur at each 
magnetic breakdown junction. If a quasiparticle crosses 
N magnetic breakdown junctions, one has to consider 
2 N possible trajectories with their correct probabilistic 
weightings, and this complicates a direct computation of 
magnetoresistance since one has to sum over trajectories 
with arbitrarily long path lengths and hence arbitarily 
large values of N. 

In order to calculate how magnetic breakdown can 
change the connectivity of the orbits^, and thus the resis- 
tivity, Falicov and Sievert in 1965 provided a models that 
generalised the Chambers path integral method 8 to in- 
clude magnetic breakdown at a finite number of points. 
This was used to evaluate the magnetoresistance in a 
number of model Fermi surfaces, although the calcula- 
tions were performed numerically. 

More recently, the problem of the effect of magnetic 
breakdown on the magnetoresistance in a network of cou- 
pled orbits was treated analytically^, though the focus of 
this work was on interlayer transport. This is because 
the newer calculations are motivated by recent experi- 
ments on quasi-two-dimensional organic metals^— in 
which magnetotransport experiments have been particu- 
larly illuminating 14 . In this paper we use the techniques 
developed in Ref. HI (see also Ref. [To] ) to derive an analyti- 
cal expression for the in-plane components of the resistiv- 
ity tensor in the Falicov-Sievert model. This expression 
allows the underlying physics of the magnetotransport to 
be more clearly extracted and will allow a more conve- 
nient comparison with experimental data. 



II. INTRODUCTION TO THE MODEL 

In this paper, we consider the quasi-two-dimensional 
Fermi surface shown in Fig. HJa) in which the dispersion 



is given by E{k) 
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FIG. 1: (a) The Fermi surface in the k x -k y plane showing 
the points where magnetic breakdown can occur which are 
at (k x , k y ) — (±fcp sin£, ±fcp cos£) (these are called magnetic 
breakdown junctions). The azimuthal coordinate of a quasi- 
particle at the point labelled QP is cp. (b) An open orbit 
(comprising the repeated traversal of the B\ section), (c) 
Closed orbit (the a-orbit, comprising the repeated traversal 
of A\ and ^b). (d) Breakdown orbit (the /3-orbit, comprising 
Ai-Si-A 2 -S 2 )- 



and g(E) = m*/d±?rh is the density of states. The 
quasiparticle orbits lie in a plane perpendicular to the 
magnetic field B and orbit the Fermi surface with angular 
frequency lo c given by lo c — eB/m* if B is perpendicular 
to the quasi-two-dimensional planes. Writing tp = ui c t + 
tpo, we have that 
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and so the conductivity tensor can be written as 
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is the effective mass, is the Fermi wave vector, d± 
is the interlayer spacing and the interlayer hopping t± 
is small (t± <C hk-p/d±). The Fermi surface consists of 
a cylinder with volume Trkp x and thus the number 
density of electrons is n — kp/2irdj_. The Brillouin zone 
edges are placed at k y = ±fcpcos£. Quasiparticles orbit 
around the Fermi surface with constant k z when B lies 
along the interlayer direction. In very low B, because of 
Bragg reflection, only open orbits [Fig. QXb)] and small 
closed orbits [the a orbit, Fig. Hfc)] occur around the 
distinct sections of the Fermi surface. In high B, mixing 
between the states on the two Fermi surface sections leads 
to magnetic breakdown at the four magnetic breakdown 
junctions. At these junctions a quasiparticle "tunnels" 
in fc-space between the Fermi surface sections, resulting 
in a single large closed orbit [the /3 orbit, Fig.QJd)]. This 
model is equivalent to the situation outlined in Fig. 1 
of Ref. lD but the notation has been chosen to link with 
the physical situation relevant to a particular family of 
organic metals&i^. 



and 



00 



e 2 m*VpT ne-f 
2i:h 2 d i 



B 



(6) 



Neglecting the small effect of the interplane warping t± 
on orbit size, and considering only the in-plane compo- 
nents of the conductivity tensor er, we proceed as follows. 
Setting v x = v-p cos (p and v y = v-p sin tp, and assuming 
full magnetic breakdown (so all the orbits just go round 
the cylinder with no Bragg scattering), the integral can 
be evaluated to give 
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and inverting this tensor yields the resistivity tensor 
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or equivalently 



III. MAGNETORESISTANCE WITH NO 
BRAGG SCATTERING 

In this section we review the standard calculation of 
magnetoresistance in the case of full magnetic breakdown 
(i.e. p = 1) which is equivalent to assuming a very high 
magnetic field and ignoring any Bragg scattering. In this 
case we consider only the effect of a single closed break- 
down orbit [i.e. the case of Fig. [ljd)]. The Boltzmann 
equation gives the conductivity tensor 
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as an average of velocity correlations over the Fermi sur- 
face, where vp is given by 
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This result is exactly what one expects from the sim- 
ple theory of metals: no magnetoresistance (because 
Pxx = Pyy = ^q 1 is not field-dependent) and a simple 
Hall effect p xy = —p yx = RhB where i?n = 1/ne. For 
later reference, we will denote the simple form of the re- 
sistivity tensor with no magnetic breakdown in Eq. (|SJ) 
as p°. 



IV. MAGNETORESISTANCE INCLUDING 
MAGNETIC BREAKDOWN AND BRAGG 
SCATTERING 

A treatment fully including both magnetic breakdown 
and Bragg scattering can be performed using the method 
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described in Ref. [9|- Setting v a = v a /vp with v x 
and v y = sin ip, Eq. (Q| becomes 
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and with the definitions of the functions E^; (x) 

E±(z) = ^(x)e ±a; /T, 
the conductivity can be written 
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What makes Eq. (fT2")l challenging to evaluate is that the 
integrand changes depending on the path taken by the 
quasiparticle which, at each magnetic breakdown junc- 
tion of the orbit, can either undergo magnetic breakdown 
tunneling (with probability p = e~ B °/ B ) or Bragg reflec- 
tion (with probability q = 1 —p)±£. The strategy for solv- 
ing this problem follows the earlier approach*^ of sep- 
arating the motion of electrons in the four different seg- 
ments of the orbit and constructing recursive equations 
which encode all the information about the behaviour at 
the magnetic breakdown junctions. In this way we can 
write the conductivity tensor as 
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where the vectors A^ and A™ 1 * are given by 
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and the matrix T which encodes all the breakdown prob- 
abilities and scattering is given by 



/ a 

6 

a 

V b 



/0p?0 

q p 

q p 

\p q 



(16) 



where a = exp (—2^/7) and b = exp (— (tt — 2^/7). 

The integrals can be performed analytically and after 
much simplification (sec Appendix) the conductivity can 
be obtained as 
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where the parameter A is given by 

4|7| ( q(l + ab) cos 2 £ 



tt(1 + 7 2 ) V 1 - &<? + a(q + b(p 2 - q 2 )) 
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and the matrix X is given by 



X = 
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7 2 sin 2 £ — cos 2 £ 
— (7 2 + 1) sin £ cos £ 



~7 — (7" + 1) sm £ cos £ 
7 2 cos 2 £ — sin 2 £ 

(19) 

Equation (fTTj) demonstrates that the conductivity is a 
sum of the familiar expression assuming no Bragg scat- 
tering [i.e. Eq. (J7)] and an additional term. This addi- 
tional term is included by an amount controlled by A, 
and it is therefore useful to plot the field dependence of 
this quantity. This is shown in Fig. [21a) for the case 
in which £ = 7r/6 (chosen to match Fig. 1 of Ref. @) 
for a range of different breakdown fields (parametrized 
by 70 = luqt = eBor/m*) and plotted as a function of 
magnetic field (parametrized by 7 = uj c t = eBr/m*). 




FIG. 2: (Color online.) The parameters (a) A [Eq. lfl8]l]. (b) 
A [Eq. ([21])] and (c) II [Eq. (Ei^ ]. 

For low field A « 0:7/(1 + 7 2 ) where a — (4 cos 2 
and is independent of breakdown field 70. For interme- 
diate fields, and for 70 S> 1, A rises to a plateau given 
by A w 4cos 2 £/[7r(7r — 2£)] = a/ (it — 2£). In our exam- 
ple in which £ = 7r/6, the plateau is at A « 9/27T 2 [see 
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Fig. Ufa)]. For large fields, A decreases again and fol- 
lows A « a7o/7 2 . Thus we expect that at very low and 
very high fields the approximation of no Bragg scatter- 
ing (the single breakdown orbit) in Eq. ([7]) will work well. 
It works at very low fields because electrons travel only 
a very short distance around the Fermi surface before 
scattering and so orbit connectivity is largely irrelevant. 
It works at very high fields because then the magnetic 
breakdown probability is close to unity. At intermediate 
fields A is significant and the effect of magnetic break- 
down and Bragg scattering is important in determining 
the magnetotransport. 
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FIG. 3: The components of the resistivity tensor (plotted 
in units of ovT 1 ) calculated according to Eqs. (I20H22|) for the 



case of £ = 7r/6 and ojqt = 30. 

We now invert the conductivity tensor in Eq. (|I7j) to 
obtain an expression for the resistivity tensor, and find 



II + Acos2£ -7 + A sin 2$ 
7 + A sin 2£ II - A cos 2£ 



where 
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Equations (|20l - |22l) are the main results of the paper. The 
quantities A and II are plotted in Fig. [U[b) and (c). In 
low field (7 <C 1) A » 07 and II w 1. In large field 
(7 > 1) A « a jo and II w 1 + cryo- These results can be 
used to show that in low field the resistivity tensor is 



p = p° + er l 'yaY, 

while in high fields the resistivity tensor is 
p = p° + cr ( 7 1 7ooY, 
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where the matrix Y contains only geometric factors and 
is given by 
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. DISCUSSION 
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The components of the resistivity tensor obtained us- 
ing equations (f20W22"]l are shown in Fig. [3] for the case of 
£ = 7r/6 and uj t = 30. At high magnetic field the diag- 
onal components of the resistivity tensor saturate (to a 
value linear in 70), while the off-diagonal terms approach 
the asymptotic values proportional to ±7 + 7oct sin 2£ [in 
agreement with Eq. As the magnetic field is de- 

creased all resistivity tensor elements go through a max- 
imum that becomes sharper with increasing breakdown 
field 70. At lower fields they decrease linearly. The posi- 
tion of the maximum is controlled predominantly by the 
behavior of A plotted in Fig. [2j 

The components of the resistivity tensor are plotted 
again in Fig. [H though this time for a range of values of 
breakdown fields (parametrized by 70 = wot). [The plot 
of p xx is identical to that of Fig. 1(c) in Ref. [3, demon- 
strating the agreement of the analytical expressions in 
Eqs. (|2"DH2"2"|) and the earlier numerical work.] The differ- 
ence between p xx and p yy reflects the fact that for our 
considered geometry in Fig. [4] (£ = 7r/6) the open orbits 
which occur due to Bragg scattering are efficient at car- 
rying current in the y-direction. This difference vanishes 
when £ = 7r/4. 

It is worthwhile to show that the off-diagonal compo- 
nents of the resistivity tensor do indeed obey Onsager 
symmetry, i.e. p xy {B) = p yx {—B). Onsager symmetry 
is dependent upon perfect microscopic reversibility (with 
the appropriate sign change in the magnetic field)^, and 
we find that this is obeyed even with the presence of mag- 
netic breakdown junctions, despite the apparent "ran- 
domization" of the electron trajectory which occurs at 
magnetic breakdown junctions. Reversal of the magnetic 
field changes the way in which the sections of the Fermi 
surface are connected at the magnetic breakdown junc- 
tions and the result of this is that the connectivity matrix 
r transforms into its transpose. This has the effect of 
making A invariant under a sign change of the magnetic 
field [and is the origin of the modulus sign of the factor 
of 7 in the numerator of Eq. (|I8I) ]. The net consequence 
of this is that Onsager symmetry is preserved. 

In summary, in this paper we have derived an ana- 
lytical expression for the in-plane components of the re- 
sistivity tensor that includes the effect of all the orbits 
induced by a combination of Bragg scattering and mag- 
netic breakdown. Eq. (JTTJ) demonstrates that approach 
shows that the conductivity tensor can be written as a 
sum of a matrix representing the effect of total magnetic 
breakdown and one representing a combination of com- 
plex electronic trajectories. Our main result in Eqs. (|20f - 
I22p provides a compact form for the in-plane components 




of the resistivity tensor for this problem which are in a 
convenient form to compare with experimental data. 
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E^(x] can be simplified because v x = cos tp, v y = siny, 
so that 



E^(yj) = cos</?e 



±<ph = ^ 



K (e^ 1 
(e^ 1 



The vector in Eq. fHl) can then be written 
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Appendix 

The expression for conductivity from Eq. ()13[) is 
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and as in Ref. H the matrix r(J — T) 1 is given by 
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where r = 1 — bq, s = q + b(p 2 — q 2 ), N = r 2 
t = a 2 s 2 and w — b(qr + a 2 (p 2 — q 2 )s). The function 
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with z ± = (z ± 1 /7) 1 . The two first and two last entries 
of each vector differ by a minus sign because these parts 
of the orbit are exactly ir apart. With this identity we 
can readily calculate: 
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The initial contribution can be simplified as follows: 
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Integration of the second term then leads to 
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The final result in Eq. (|T7|) is obtained by substituting expressions for A^ ± and A„ ± from Eq. (|28|) and for A 1 ; and 
\ 2 ai from Eq. (|2"TJ)) and simplifying. 
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